indv=find(p1(1,:));

bi=[0,1;3,6;2,3;2,6;2,9;2,12;3,4;6,7;4,5;7,8;9,12;9,10;12,13;10,11;13,14];
b=bi+ones(15,2);

c=zeros(1,15);
  for j=1:15
     c(j)=(b(j,1)-1)*15+b(j,2);
     
  end